Finite-time rotation number: a fast indicator for chaotic dynamical structures 
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Lagrangian coherent structures are efFective barriers, sticky regions, that separate phase space 
regions of different dynamical behavior. The usual way to detect such structures is via finite-time 
Lyapunov exponents. We show that similar results can be obtained for single-frequency systems 
from finite-time rotation numbers, which are much faster to compute. We illustrate our claim by 
considering examples of continuous and discrete-time dynamical systems of physical interest. 
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It is noteworthy that chaos is observed in physically 
interesting systems, both in Nature and in mathematical 
models. Indeed, chaotic dynamics is now commonplace in 
such diverse disciplines as celestial mechanics [l[ , atomic 
physics [3], fluid mechanics [1], plasma physics etc. 
The main means of determining chaos is with Lyapunov 
exponents, which can be notoriously difficult to compute 
in fast and reliable ways, particularly if only experimen- 
tal data are available and the governing equations are 
unknown. Hence, an active area of research is the search 
for fast indicators (FIs), i.e., coinputational diagnostics 
that characterize chaos quickly [5|. Here we introduce 
a new FI, the finite time rotation number (FTRN), for 
determining effective transport barriers or sticky regions. 

Since Lyapunov exponents measure the exponential 
rate of divergence of nearby trajectories, FIs are usually 
based on them, examples being the finite-time Lyapunov 
exponent (FTLE) Q, the smaller (SALI) and general- 
ized (GALI) alignment indices Q, and the mean expo- 
nential growth of nearby orbits (MEGNO) [1]. These FIs 
not only determine the strength of chaos, but can detect 
invariant tori and other issues related to integrability. 
Alternatively, there are diagnostics based on frequency 
decomposition, which are well-suited for weakly chaotic 
motion with any number of degrees of freedom Q . 

A relatively new application of FIs is to quantita- 
tively characterize chaotic transport, an important basic 
physical process occurring in many contexts, from fluid 
dynamics 0, to fusion plasma confinement A 
now popular method for determining chaotic transport 
is based on Lagrangian Coherent Structures (LCS) jlli] . 
LCS represent effective barriers that separate regions of 
different dynamical behavior. This FI has been used to 
investigate time-dependent flows occurring in many ap- 
plications of fluid mechanics, such as transport in ocean 
currents [12], flow over an airfoil fli|. and chaotic mixing 
in forced tanks [14], and it has recently been introduced 
into plasma physics and used to describe plasmas turbu- 
lence [l5| and magnetic reconnection [l6| . 

In practice LCS are unveiled via computation of 
FTLEs for many points, i.e. the FTLE field. One then 
finds ridges in this FTLE field that inhibit transport be- 
tween different regions. Such ridges, although not bona 



fide invariant sets, provide effective transport barriers. 

The detection of LCS by FTLEs is versatile ~ it can 
be applied to periodic, quasiperiodic, or broadband vec- 
tor fields defining the flow. However, in two-dimensions 
when the vector field is area-preserving and contains a 
single temporal frequency, the FTRN proposed in this 
article is a better alternative. Although for tractabil- 
ity many early studies of mixing in fluids Q consid- 
ered single- frequency vector fields, physical velocity fields 
contain many frequencies or broadband turbulent spec- 
tra. Moreover, area preservation in fluids arises from the 
solenoidal approximation of the velocity field, a common 
assumption, of e.g. geophysical interest (e.g. [13, [isj). 
with varying degrees of validity. However, for analyzing 
snapshots of magnetic fields in toroidal plasma devices 
or, more generally, astrophysical or other magnetic fields 
in the vicinity of any stable closed field line 2l|, the 
time-like coordinate is a toroidal angle [l^ and the sys- 
tem contains exactly only a single frequency. Also, area 
preservation is an exact consequence of V • B = 0. Thus 
for such systems, the FTRN is a natural and better suited 
fast indicator for characterizing chaos and transport. 

For single frequency systems, we show results on LCS 
using the FTLE can be obtained by the simpler and com- 
putationally faster method based on the FTRN. As for 
the FTLE, LCS are ridges of the FTRN field computed 
from a grid of initial conditions. But, the rotation num- 
ber does not require the evaluation of spatial derivatives; 
thus a fine mesh, although desirable, is not essential and 
this substantially reduces the computational time to ob- 
tain LCS with good resolution. 

A time-T periodic dynamical system with annular 
phase space V is determined by a Poincare map, M : V ^ 
T), where each period-T is represented by one iteration of 
M . If X G 2? lies on an invariant circle , then M maps 
to itself. The rotation number (e.g. [17|) for an or- 
bit starting at Xq is = lim„_>oo H • (M"(xo) — xo)/n 
which is lifted to M and H is a suitable angular pro- 
jection. Under mild conditions on M, this limit exists 
for every initial condition Xg G and does not depend 
on xq. Consider a simple example, the rigid rotation 
Af(x) = (x + w,y), where x — {x,y) and x 6 S^. Here 
the rotation number lo ^ w. If is a rational number 
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p/q, the trajectory is a period-g orbit of the map M, 
whereas if w is irrational, then the ensuing (quasiperi- 
odic) orbit covers densely the circle . The FTRN is 
the time-iVr truncation for the flow and corresponds in 
the preceding definition to N iterations of M . Below wc 
use Wiv(xo) := H • (M^(xo) - xo)/iV, with only N 
iterations. In general, cjjv, like any truncation, depends 
on the initial condition. While the infinite-time rotation 
number is not defined for chaotic orbits, which do not lie 
on any 5^, the finite-time counterpart exists for any or- 
bit. Roughly speaking, wat measures the average rotation 
angle swept out by a trajectory over a time interval NT, 
and thus conveys information about the local behavior of 
trajectories, just as the FTLE does (which measures local 
rates of contraction or expansion) . We identify LCS with 
ridges of the FTRN. Note, minor ridges of the FTRN are 
finite pieces of invariant tori, which appear because the 
full rotation number has not been calculated. 

In order to illustrate how the FTRN detects LCS, we 
consider three examples and compare each with FTLE 
calculations. Essentially the same results are obtained 
with FTRN, but with lower computational cost. One 
example is a time-periodic two-dimensional fluid flow, the 
second is a discrete-time map of a flow used for passive 
advection, and the third is a magnetic field-line map. In 
all cases a uniform grid of 800 x 800 points is advanced. 

1. Periodic double gyre flow. We consider a two- 
dimensional fiuid fiow with a stream function, 

'4){x, y)= A sin[7r/(a;, t)] sin(7r?/), (I) 

where /(a;, = a{t)x'^ + h{t)x , a{t) = e sin(27rt/T), h{t) ^ 
I - 2esin(27rt/T), 2? := {0 < a; < 2, < j/ < 1}, and A 
is the maximum horizontal velocity, u [l9j . The velocity 
field V = (w, v) is given hy u = —dip/dy, v = d^/dx. 

For e = the flow is integrable, with equilibrium points 
A : (1/2,1/2), B : (3/2,1/2), C, : {,x,i,y,i), z = 1,...6, 
where Xd £ {0, 1, 2} and yd G {0, 1}. The points A and 
B are centers, whereas Ci are saddles connected by het- 
eroclinic trajectories. The latter are boundaries of two 
gyres surrounding A and B, with clockwise and counter- 
clockwise rotations, respectively. The heteroclinic tra- 
jectory 'H connecting C3 : (0, 1) and C4 : (1, 1) separates 
two distinct gyres and thus is a natural place to focus 
attention when looking for LCS. 

For e 7^ the flow is time-dependent and noninte- 
grable, yet T) remains invariant. As is well-known, the 
former heteroclinic connections are structurally unsta- 
ble; upon perturbation an entanglement of stable and 
unstable manifolds with concomitant horseshoe dynam- 
ics appears. The vertical line % is no longer invariant, 
but can be thought of as roughly separating two gyres 
with time- varying amplitudes: % oscillates in the hor- 
izontal direction with amplitude ~ e (for small e) and 
frequency 27r/T. The T-periodicity permits a time-4r ro- 
tation number, a;4(x) = X^j^^ ^(xi)/4, where tan0(x) = 
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FIG. 1: (color online) (a) Time-4r Lyapunov exponent and 
(b) time-4r rotation number for the double gyre system, with 
period T = 10, amplitude A — 0.1 and forcing strength e = 
0.25. (c) and (d) depict the Lagrangian Coherent Structures 
corresponding to ridges of (a) and (b), respectively. 

(xA.B — x)/{yA,B — 2/) is a rotation angle around gyre A 
or B, for the time-4T stroboscopic map of this noninte- 
grable system. This defines the projection 11. Similarly, 
we can compute the corresponding time-4r Lyapunov 
exponent, which requires computing five nearby orbits, 
necessary to evaluate spatial derivatives, instead of only 
one orbit for the FTRN. For this reason, computation of 
FTRNs is at least five times faster than FTLEs. 

The results of the FTLEs and the FTRNs for the dou- 
ble gyre system are depicted in Figs. [TJa) and (b), re- 
spectively, where blue (red) depict low (high) values of 
the corresponding quantity. The ridges of both, i.e. the 
crests of higher values are shown in Figs. [TJc) and (d) 
for FTLE and FTRN, respectively. The pictures are in- 
deed very similar, notwithstanding the wide difference 
in the CPU-time necessary to produce them. Moreover, 
Figs. [TJc) and (d) reveal the existence of LCS for the 
nonintegrable system. For small values of y, the LCS ap- 
proach the oscillating vertical line H that separates the 
gyres. In fact, despite the absence of well-defined sta- 
ble and unstable manifolds of equilibria for time-periodic 
fiows, the LCS displayed by Fig. [1] are quasi-invariant: if 
a passive scalar (tracer) were put on such LCS, it would 
be advected by flow and remain in the vicinity of the LCS 
for a long time (on the order of the experiment duration). 
In other words, even though the LCS no longer separate 
the gyres for an arbitrarily long time (there may be a 
small transverse flux), trajectories starting on the left- 
handside (righthandside) chiefly remain in the lefthand- 
side (righthandside). In practical terms, however, this 
suffices to characterize an effective transport barrier. 
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FIG. 2: (color online) (a) Time-IOOT Lyapunov exponent and 
(b) time-lOOr rotation number for the advection map ((2]) with 
mixing parameter a — 0.25. (c) and (d) depict the Lagrangian 
Coherent Structures corresponding to ridges of (a) and (b), 
respectively. 



2. Advection of a passive scalar. In Example 1., in- 
vestigation of LCS in a time-T periodic flow required 
the numerical evaluation of the time-T stroboscopic map 
M . Sometimes when researchers consider advection of a 
passive scalar by two-dimensional flows, M is given by 
replacing the mixing action of a flow by an explicit area- 
preserving map. An example ,20;] of this is 



Vn+l = Un 



- a sin(27rj/„) (modi), 
asin(27ra;„) (modi) 



(2) 



The map ([2]) corresponds to a velocity field that is the 
superposition of two sinusoidal shear flows in the x and 
y directions. In this case, the flow shear reverses sign 
along some shear less curve, as is the case for zonal flows 
of geophysical, atmospheric, and plasma physical interest 
[Tct |2l|. The map (O is symplectic and represents a 
Hamiltonian system for any value of a. The fixed points 
in the torus T> = [0, 1) x [0, 1) are the centers P : (0, 1/2), 
Q : (1/2,0) and the saddles R : (0,0) and S : (1/2,1/2). 

Since the system is nonintegrable for a ^ 0, the stable 
and unstable manifolds stemming from the saddle points 
R and S intersect in a heteroclinic tangle and there are 
chaotic orbits that do not line on continuous invariant 
circles. This structure is responsible for the mixing effect 
of the chaotic advection. Nevertheless, this chaotic layer 
acts as an effective transport barrier separating the two 
gyres with invariant tori encircling P and Q. 

This is clearly seen after computing the time-lOOT ro- 
tation number for M (with 11 projecting onto the x-axis) 
and then extracting the corresponding ridges with high 
values of ujioo- These ridges trace out quasi-invariant 
sets that shadow the heteroclinic connections, especially 



in the vicinity of the saddle points R and S, reinforcing 
their interpretation as LCS. The computation of FTRNs 
is particularly fast for maps, so calculation of wioo is not 
difficult, but the results are essentially identical for ujiq. 

In Figs. (Ha) and (b) we depict the FTLE and FTRN, 
respectively, for orbits of the advection map ^ for 
a = 0.25, whose ridges (points with largest relative val- 
ues) are shown in Figs.[2l^c) and (d), respectively. Both 
diagnostics indicate that the LCS are quasi- invariant sets 
about a chaotic separatrix layer. This layer acts as a 
transport barrier that separates quasiperiodic curves en- 
circling the centers P and Q; the layer being a ridge 
implies these are LCS of this system. 

3. Magnetic field line map. Magnetic field lines in the 
equilibrium states of toroidal magnetic plasma confine- 
ment devices, such as tokamaks and stellerators, are or- 
bits of a one degree-of-freedom integrable Hamiltonian 
system, where a toroidal-like angle plays the role of time 
(e.g. 21|). In the simplest case, canonically conjugate 
variables are the spatial coordinates of an annulus and 
points (xmUn) represent the n-th field line intersection 
with a surface-of-section at a fixed value of the toroidal 
angle. A similar situation arises in the vicinity of any 
stable (elliptic) closed magnetic field line. 

Perturbations due to external electric currents or in- 
ternal MHD instabilities break symmetry and render the 
system nonintegrable, giving rise to chaotic field lines. 
Since the magnetic field configuration here is strictly 
static in time, the term chaos means that two infinites- 
imally close field lines Lyapunov exponentiate as they 
wind around the torus. 

An example of a nonintegrable system that models 
such field lines is provided by the so-called tokamap [l^] : 



Vn+l = Vn- T^Z 
= X 



L_ yn+i sin(27rxn) 
2tt 1 + yn+i 

1 L cos(27ra;„) 



q{yn+i) 27r (1 + y„+i)2' 



(3) 



where L is a parameter measuring nonintegrability (pro- 
portional to the perturbation strength) and 1/g = (2 — 
y){2 — 2y + 2/^)/4 is the inverse of the rotational trans- 
form of the field lines. The value of yn+i is obtained by 
applying Newton's method to ^ at each map iteration. 

For i = 3/27r the FTLE and FTRN (with H again 
the x-projection) of the tokamap are shown in Fig. [Hlja) 
and (b), respectively, with (c) and (d) depicting the cor- 
responding ridges. In both cases the latter correspond, 
as in the advection map, to the homoclinic intersection 
of manifolds of the saddle fixed points of the tokamap 
([3]), defining a thin chaotic layer originating from such 
intersections. The separatrix layer is again a LCS, which 
indicates a transport channel along the layer, yet restrict- 
ing diffusion across the barrier. 

In conclusion, here we have proposed the FTRN as a 
FI, which for physical systems with a single period is su- 
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FIG. 3: (color online) (a) Time-IST Lyapunov exponent and 
(b) time-lST rotation number for the tokamap ((3| with L = 
4. (c) and (d) depict the Lagrangian Coherent Structures 
corresponding to ridges of (a) and (b), respectively. 



perior, being faster and simpler, to the FTLE. The three 
examples treated demonstrate this point. We note, how- 
ever, that the speed of calculation of both FTRNs and 
FTLEs is achieved at the cost of sacrificing the detailed 
picture of transport provided by an analysis of turnstiles 
through cantori 



e.g. [23|), which requires a search for 



periodic orbits [17|. In future studies we propose relax- 



ing the single frequency limitation by a more detailed 
frequency analyses (e.g. and considering statistical 
analyses of the FTRN akin to that of FTLE (e.g. [s]). 
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